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A.  A.  SAMARSKIY 


1 .  The  broadening  of  the  scope  of  difference  methods  and 
the  associated  increase  of  demands  on  them  stimulate  interest  in 
formulating  general  principles  of  obtaining  difference  schemas  of 
a  definite  quality  for  problems  of  mathematical  physics.  The 
possibility  of  assigning  to  the  same  differential  equation  an 
infinite  set  of  difference  schemas  defined  on  the  same  template 
and  having  the  same  order  of  approximation  makes  the  problem  of 
choosing  the  desired  schema  a  very  complex  one.  It  is  required 
of  any  numerical  method  that  it  supply  an  approximate  solution  of 
the  problem  with  preset  accuracy  ^  >  0  in  a  finite  number  of 
operations.  Moreover,  it  is  required  that  the  schema  be 
universal  (suitable  for  a  rather  broad  class  of  problems), 
homogeneous,  stable,  and  economical  (more  precisely,  that  the 


numerical  algorithm  used  to  solve  the  difference  equations  be 


economical ) . 


-  ^  ^  ’  *■  'r  •>  * 


]/  n  l a 


All  these  requirements  vie  with  each  other. 

Homogeneity  of  the  schema  Cl]  means  that  the  difference 
operator  is  defined  by  the  same  formula  at  all  grid  nodes  for  any 
coefficients  and  right  side  of  the  equation  of  a  given  functional 
class,  and  also  for  an  arbitrary  grid. 

Below  we  examine  only  homogeneous  schemas. 

2.  Approximation  error  is  used  as  one  of  the  a  priori 
characteristics  of  a  schema.  The  greatest  possible  order  of 
approximation  of  a  solution  is  determined  under  the  assumption  of 
adequate  smoothness  of  the  solution.  A  local  Cat  a  point) 
approximation  is  not  generally  a  very  good  characteristic  of  a 
schema.  The  negative  norm  Cor  “integral  norm” >  Cof  the  form 
||  til  A-l  ■  Vc A“1  <f  ~  Y~ )  ,  a  =  A»  positive ' definite  in  case  of  the 
difference  equation  Ay  ■  H'  ) ,  is  a  natural  norm  for  estimating  the 
error  of  approximation. 

In  the  case  of  linear  stable  schemas,  the  error  of  the 
schema  depends  continuously  on  the  error  or  approximation,  so 
that  the  error  of  approximation  determines  the  error  of  the 
schema.  For  nonlinear  equations  Ce.g.,  for  gas  dynamics 
equations)  such  estimates  are  not  available  Cso  far  there  is  no 
proof  of  stability  of  any  of  the  difference  schemas  for  gas 
dynamics) ;  so  here  the  criterion  of  approximation  plays  more  or 


2 


leas  a  formal  role  and  may  give  an  improper  impression  aa  to  the 


accuracy  of  a  schema  owing  to  strong  discontinuities. 

Note  that  comparison  of  schemas  in  terms  of  order  of 
accuracy  is  meaningful  only  for  sufficiently  small  steps  h  and 
'T' ;  in  this  connection  the  higher  extent  of  step  is  multiplied  by 
the  maximum  of  the  derivative  of  a  solution  of  higher  order.  But 
in  practice  very  coarse  grids  are  used  where  asymptotic  estimates 
may  not  work.  It  can  happen  that  a  schema  of  first-order 
accuracy  in  an  actual  grid  is  more  accurate  than  a  schema  of 
second-order  accuracy. 

Thus,  the  criterion  for  choosing  schemas  should  be  stated 
as:  the  schema  should  give  sufficient  accuracy  on  actual  grids 

for  the  class  of  problems  under  consideration. 

From  this  point  of  view  one  should  take  into  account 
qualitative  considerations  as  well  as  mathematical  ones  (if  there 
are  any ) . 


3.  Conservativeness  of  the  difference  schema  is  one  such 
requirement . 

What  i3  conservativeness? 

For  definiteness,  we  will  deal  with  problems  of  mechanics  of 
a  continuous  medium. 
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In  setting  about  solving  such  problems  in  some  region  by  the 
method  of  finite  differences,  one  introduces  a  grid  in  the  region 
G  and  replaces  the  differential  equations  by  difference  equations 
for  the  grid  functions.  As  a  result  one  gets  a  mathematical 
description  of  a  discrete  model  of  the  medium.  Obviously,  the 
discrete  model  should  reflect  the  main  features  of  the  continuous 
medium.  The  properties  of  a  continuous  medium  are  determined  by 
the  integral  laws  of  conservation  (of  momentum,  mass,  energy, 
etc.)  for  any  subregion  G'  in  G. 

The  differential  equations  are  corollaries  of  the  integral 
laws  of  observation.  It  is  natural  to  require  that  the 
difference  schemas  express  the  laws  of  conservation  on  the  grid. 
The  laws  of  conservation  for  the  entire  grid  fegion  should  be 
algebraic  corollaries  of  the  difference  equations.  Schemas 
possessing  these  properties  are  called  conservative. 

15-20  years  ago  the  question  of  whether  the  schema  should  be 
conservative  ("divergent")  might  have  been  a  matter  of 
discussion,  but  at  present  there  is  general  agreement  that  it 
should  be. 

Along  with  conservative  difference  schemas  one  also  uses 
conservative  differential -difference  schemas. 

The  general  method  of  integral  ratios  for  obtaining 
conservative  differential -difference  schemas  was  suggested  by 


A. A.  Dorodnitsyn  C23  and  applied  for  solving  multivariate  gas 
dynamics  problems.  Later  this  method  was  further  developed  by 
his  students  C3] . 

To  get  conservative  difference  schemas  one  can  use  the 
integro-interpolational  method  or  the  method  of  balance  [1] . 

A  difference  operator  on  the  space  of  grid  functions  should 
preserve  the  basic  properties  of  the  differential  operator 
defined  on  the  space  of  functions  of  continuous  argument.  Such 
properties  in  the  linear  case  are  self -conjugacy  and  (sign-) 
definiteness  of  the  operator.  Below,  with  the  example  of 
Laplace-operator  schemas  in  an  arbitrary  regions,  we  prove  a 
relationship  between  self -conjugacy  and  conservativeness.  A 
nonconservative  operator  is  also  nonself -conjugate . 

The  disadvantages  of  nonconservative  schemas  are  not 
removable.  Condensation  of  the  grid,  which  one  ought  to  be  able 
to  count  on  in  order  to  increase  accuracy,  may,  in  the  case  of  a 
nonconservative  schema,  even  increase  the  error  of  the  schema. 
Example  in  Cl]  confirms  this  for  the  problem  <ku')'  *  0, 

0  <  x  <  1,  u(0)  *  1,  u(l>  =  0. 

The  nonconservative  schema 

( bi  (y,+ 1  ~  l/i)  —  a Au,  —  Ui-Jlh2  =  0,  1  =  1.  2 . .V  —  J , 


where  =  ki  -O,  25  (ki  +  i  -  kj.-i>,  bi  *  ki  ♦  0,25  <ki  +  i  -  ki+i), 
diverges  in  case  of  piecewise  constant  coefficient  k(x)  =  kj  for 
x  <  ^  and  k(x)  =  Kji  for  x  >  ^  ,  and  this  divergence  is  quite 
peculiar:  the  solution  <xi)  of  the  difference  problem  has,  as 

h-*0,  a  limit  u(x)  not  equal  to  the  exact  solution  u(x>  of  the 
differential  problem. 


The  limit  function  u(x>  is  a  solution  of  the  problem  with 
the  additional  condition:  at  the  point  x  =  ^  of  discontinuity  of 
k(x)  there  is  put  a  source  of  power  depending  on^  ,  kT,  kjj, 

equal  to  zero  only  for  ki  *  ku  and  becoming  infinity  for 
certain  ki,  >'11,1^  .  Thus,  verification  of  the  accuracy  by 

condensing  the  grid  may,  in  this  case,  lead  to  the  improper 
conc_usion  that  the  nonconservati ve  schema  converges. 


Note  that  any  difference  schema  generates  fictitious  sources 
(drains).  Indeed,  let-'A-y  * ‘f  =  0  be  a  difference  schema,  u  the 
exact  solution  of  the  differential  equation  Lu  *  f  -  0.  The 
discrepancy-/]-  u  +  Cf  =  'f'  is  the  error  of  approximation;  it  may  be 
viewed  as  the  density  of  the  fictitious  sources. 


If  the  schema  is  conservative,  thenl^(x)  is  an  oscillating 
function,  so  that  always  we  have 


(y,  1)  =  ^  i  (r)  h  -►  0  hpn  h  -►  0. 


Indeed,  for  <ku')'  =  -f<x>,  0  <  x  <  1,  the  balance  equation 
has  the  fora 

i  i 

k u  -f  j  f(x)dx  =  0. 
o  o 

The  corresponding  conservative  schema  (a(x)yx>x  = 
x  =  ih,  possesses  the  analogous  property  on  the  grid: 

if  v  —  wi  +  (=•  1)  =  0,  w%  =  a^X'  i  =  at  y*  . 

For  =  -A-u  *  we  get 

(y,  l)  =  (i2U1).v—  (aUilj-f-ls,  l)  =  0(/r), 

if  the  schema  has  second  order  of  approximation. 

4.  The  conservativeness  of  the  homogeneous  difference 
schema  for  the  equation 

div  (k  grad  u)  =  — /  (x),  x  =  (x1 . xf) 

is  necessary  and  sufficient  for  convergence  in  the  class  of 
discontinuous  coefficients.  For  the  one-dimensional  case  this 
assertion  is  proved  in  [1].  In  the  multivariate  case,  for  an 
arbitrary  region  the  convergence  of  the  conservative  schema  with 
rate  0(V~h>  was  proved  in  [43.  The  conservative  schema  will  be 
given  below. 


-V* jf«JT  .  ■ 


Consider  the  Dinchlet  problem  for  the  Poisson  equation  in 


the  reg 


ion  G  *  [^ 


on  the  plane  x  =  fxi ,  X2> : 


—  '  —  — 
dx\  oz\ 


— f(x),  x£G.  a  jr  =  ate). 


For  this  equation  the  balance  equation  is  satisfied: 


5  ^ds+\  f  (*)<**= o. 

r  a 

Let  Cor*  =  oJh  ♦  y  h  ^  the  grid  in  G  ♦  F 


described  in  C51 .  It  is 


uniform  (with  steps  hi  and  xi  and  h2  in  X2>  at  strictly  interior 
nodes  x  £  ob  and  nonuniform  at  boundary  nodes  x  £  «  . 


First  let  us  consider  the  well-known  five-point  schema  C5, 


Ay  =  yA*,  -f  y^,  =  —j  (x) 


(&.) 


at  regular  nodes. 


Ay  =  A’y  +  Aly  =  ytii)  -j-  yt4t  =  —f  ( x ) 


at  nonreguiar  nodes,  where 


M'-rv.-tf »,= 0,5 


and  h*  is  the  distance  of  the  nonreguiar  node  x  6  o<~>  ««  from  th« 
boundary  node  x(+1i  -f  x<-l|  ^  h  13  a  boundary  node 

and  x("tei  >  is  an  interior  node,  then 


V  v  •/«••*»*  •  ,  *  .  • ,  .  -  -  s’  .*  . 


*.  A 


We  will  show  that  this  schema  is  not  conservative. 


Introdu 


cing  the  notation 

(y,  v)m=Zs (x) ^(x) A-A  +  -.yw n^'  (s') 

w©  get. 

(Ay,  1).  =  2  (y*fcynh i)  +  A, 

A  =  2  l\yfh  (A.  -  A.)  +  A„y/t;  (ft,  -  A)], 

*tr** 

where  yn^  =  y*^,  if  x<+io<)  is  a  boundary  node,  y ^  =  -v- 

is  a  boundary  node,  c*  =1,2. 

Comparing  (2)  and  <3),  we  see  that 

—  (y«A>  4-  y»  hx)  -r  (r-  1).=  — A,  A=^0, 

Xfce/* 

i.e.,  the  schema  Cl)  is  not  conservative. 

We  find  similarly  that  _yl-  is  nonself -con j ugate,  i.e., 

(Ay,  v).^(y,  Ay)., 

where  y  and  v  are  arbitrary  functions  defined  on  ^  and 

vanishing  on  the  boundary  of 

One  can  find  a  nonconvex  region  and  a  grid  such  that  —  _y\_ 
is  not  positive  definite. 


Using  the  integro-interpolational  method  Cl],  one  can  get 


the  five-point  conservative  schema  C4 J : 

A y  =  —f(x),  x£*k,  y  !»*?(*)•  (&>) 


At  regular  nodes,  _/L  has  the  usual  form  (2),  while  at 
nonregular  nodes  it  is  given  by: 


Ky  =  -X’y  -4*  A?y> 


<f7j> 


where 


f  »t+,,)  - 

y  v  —  vl~h)\. 

A 

i>i  ) 

.  _  i  | 

fy(”1,!-V  y- 

-*y  ~ 

V  **, 

(j/x,  —  y*,)  =  y*>*r 


CD 


if  x^^-t  )  and  x**1!*  are  boundary  nodes.  It  is  not  hard  to  see 
that 

A'=xA'-  ■ 

where  X.  and  X  have  the  form  (4) . 
i  x 

On  the  set  _A_  of  grid  functions  yCxi,  X2>  defined  on  the 
grid  ^  and  equal  to  zero  on  the  boundary  of  , 

introduce  the  scalar  product 


(y.  v)=  2 

Ifeui* 


The  above-defined  operator  _/L_  is  self -conjugate  on  _Al 


while  JL  is  positive  definite  on  an  arbitrary  grid  in  the  caae 
of  an  arbitrary  region. 


At  nonregular  nodes,  JL  has  zero  order  of  approximation: 


Au  —  lu  =  0(l)  for  x£<»“ 


However,  the  schema  (6)  has  second-order  accuracy 

| y  — B|c«0(lfcl»).  | kp*fc*  +  fc|. 

For  the  equation 

div  (A  grad  u)  =  — /  (x),  !£(?,  u  |r  =  u  (x),  *(x)>c,>  0 

the  corresponding  conservative  schema  has,  at  regular  nodes,  the 
form 

Ay  =  4-  =  — /(x). 


at  nonregular  nodes  the  form 

Ay  =  +  (fl,y1;)i:  =  —  f(x). 


where 


,(+l«)  __  | 


*; 


if  x(+l  5  ,  cK  =  l,  2.  The  coefficients  a^  are  chosen 

according  to  £5] . 
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S.  Now  let  us  turn  to  the  problem  of  conservative  schemas 
for  equations  of  magnetic  hydrodynamics  C7,  8] . 


We  first  consider  gas  dynamics  equations  which,  in  Lagrange 
variables  in  the  plane,  have  the  form: 


dp 

dr  _ 

dr,  _  dv 

(/aJ 

‘  ox  • 

~oi  V' 

~oi  ox  * 

dv 

P~a 7> 

P=P 

*). 

OI,d 

where  t  is  time,  r  is  Euler's  coordinate,  x  is  the  Lagrange  mass 
coordinate,  ^  is  specific  volume,  p  is  pressure,  £  is  intrinsic 
energy,  v  is  velocity. 


The  nondivergent  equation  <lli>  can  be  written  in  entropy 


form 


«>9> 


and,  by  using  equations  (10),  can  also  be  transformed  into  the 
divergent  form 


Olj) 


expressing  the  law  of  conservation  of  total  energy.  The  three 
systems  of  equations  (10),  (llo<  > , ^  =  1,  2,  3,  are  equivalent. 

It  is  ordinarily  assumed  (see  C9,  Ch.1113)  that  to  obtain  a 


conservative  schema  it  is  sufficient  to  approximate  the  three 


basic  laws  of  conservation  ( balance )- -of  mass,  momentum  and  total 


energy.  However,  in  this  case  the  balance  equations  may  be 
violated  for  the  individual  forms  of  energy-intrinsic,  kinetic. 
The  amount  of  disbalance  in  the  case  of  strongly  varying 
solutions  may  become  comparable  with  the  total  energy. 

We  will  call  a  difference  schema  completely  conservative  if 
both  the  laws  of  conservation  of  mass,  momentum  and  total  energy 
as  well  as  the  individual  balance  of  energy,  kinetic  and 
intrinsic,  are  valid.  To  obtain  such  schemas,  in  addition  to 
requiring  conservativeness,  a  formal  condition  is  imposed:  the 
difference  schema  should  possess  the  same  property  as  the  system 
of  differential  equations  (10),  (11c*;  ),  namely  the  "divergent" 

difference  equation  for  energy  (the  analog  of  (II3))  should  be 
transformable  both  to  the  nondivergent  form  (to  the  analog  of 
(111))  os  well  as  to  the  entropy  form  (llj;). 

A  five-point  family  of  schemas  is  examined  on  the  six-point 
template  (xi  =  ih,  t  =>  j  )»  <xi,  t;  .  ),  (Xi+.  ;tr  ), 

J  W  " *  si 

(xi *,  t^ * , ) .  By  requiring  the  indicated  conditions  to  hold, 
we  get  a  one-parameter  family  of  completely  conservative  schemas 

v.  =  —/>/',  r,  =  rtl=sv^*\  «,  =  —  p''^,  (&0 

-A 

where  we  have  used  the  notation  f  (  ^- )  =  oL  f  *  (l-=*)f,  f  =  f-, 

*  ;  a 

f  -  fJ>1,  ft  *  (f-f)/  T »  an  arbitrary  number.  Here 

v  =  V(xi,  ).  P  »  P<xi-1/2,  *  ),  n  -  ^  t-  ), 

Z  -  Z ;  (xi-i/2#  ^  ) . 
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The  nondivergent  energy  equation  £ t  -  iS 

transformed  to  the  divergent  form 


(*  ”  ~T)t  =  — <P‘*’  (— 1 )  p(-V,  =  p(x-  h), 


ct&> 


if  one  uses  the  motion  equation  Vt  =  ~Px<e,<)  and  the  formula  for 
difference  differentiation  of  the  product  (f(-l)V)x  *  fvx  +  Vfx. 


Indeed,  we  have 


=  _p"> fj*- =  -(p“’  (-1)  t,,a'5)),  +  Pf^*’  = 

=_<p-»  ^ — 1)  i;  ••*'),  -  =  — (P**‘  (-1) ~  0>5 


Whence  ( 13) . 


The  schema  (12)  has  approximation  0( T ~  h2)  for  any  o< 
Taking  *  0,5,  we  get  a  unique  schema  OC'T'2  ♦  h2>  (it  is  also 
derived  by  somewhat  different  arguments  in  [10] ) . 


7.  Now  we  turn  to  a  system  of  equations  of  magnetic 
hydrodynamics.  Let  H  =  Hx(x>  and  E  =  E^(x>  be  the  nonzero 
components  of  magnetic  and  electric  fields. 


In  Lagrange  variables  the  system  has  the  form 


where  F  is  the  Lorentz  force,  Q  is  the  Joule  heat,  6"  is  the 


electric  conductivity, 


The  energy  equation  can  be  replaced  by  the  equation  in 
entropy  fora 


d*  dr 

at  — p-zt+Q 


and  by  the  divergent  equation  for  total  energy 


d  /  .  V*  ,  8*r\  d  r  /  ,  B'-\  EH~\ 

’3r(,+  T+  — +  «r)~sr J. 


The  force  F  can  be  also  written  in  the  divergent  form 


Thus,  we  obtain  six  equivalent  systems  of  equations. 


We  will  require  that  the  analogous  equivalence  also  hold 
for  the  difference  schemas.  Examining  the  eight-parameter  family 
of  schemas  on  the  above-indicated  six-point  template,  we  get  the 
two-parameter  family  of  schemas  0(T  *  h2)  and  the  lone  schema 
0<T2  ♦  h2>  (for  =<  =  P  =  0,5): 


p,  =  -Pr  -  (4r), ,  .  r, = tk = ^-5*. 

(Srj,  =  £■/',  = 

t,  =  E'*'  -|- (5(t1)t,(+ 1)£(^1))(0-*’ 


where 


/.=  0.5 (/  +  /(-!)),  }(±\)  =  j(x-h,  t). 


The  energy  equation  can  be  transformed  to  the  divergent. 


form 


(, + 1  + *  (-1 )) + = -  [(*•■ + 


4r.  jt* 


It  should  be  noted  that  conservative  schemas 
approximating  the  total  energy  equation  may  poorly  approximate 
the  equations  for  the  intrinsic  energy  and  for  magnetic  field. 
This  is  no  less  dangerous  than  violation  of  the  law  of 
conservation  of  total  energy,  and  it  may  lead  to  an  erroneous 
computation  of  temperature,  particularly  to  such  a  nonphysical 
effect  as  the  diminution  of  the  temperature  of  some  mass  of  gas 
in  the  process  of  compression  or  in  the  presence  of  Joule 
heating.  The  disbalances  arising  here  cannot  be  eliminated  by 
spatially  condensing  the  grid. 


In  the  case  of  implicit  schemas,  to  solve  difference 
equations  one  uses  iterative  methods,  which  may  generally  violate 
the  conservativeness  of  the  schema.  Hence  the  iterations  must  be 
carried  through  for  the  prescribed  accuracy  which  characterizes 
the  amount  of  disbalance.  In  this  case,  to  simplify  the 
calculations  it  is  advisable  to  use  the  energy  equation  in 


nondivergent  form 


For  nonlinear  difference  schemas  there  is  an  amply 
developed  theory  which  allows  general  methods  to  be  formulated 
for  obtaining  set-quality  difference  schemas. 

Let  us  point  out  two  rather  general  methods: 

1)  The  method  of  regularization  of  difference  schemas 
based  on  using  the  class  of  stable  schemas  and  the  possible  of 
varying,  without  affecting  stability,  one  of  the  operators  of  the 
schema  so  as  to  satisfy  the  collateral  constraints  of  efficiency 
and  approximation: 

2)  The  method  of  summary  approximation  based  on  using  a 
new  notion  of  a  schema  and  a  new  notion  of  approximation  for 
schemas- -summary  approximation  (such  schemas  are  called 
adaptive)  . 

Both  methods  are  used  successfully,  particularly  for 
obtaining  economical  difference  schemas  in  the  case  of 
multivariate  problems  of  mathematical  physics.  A  presentation  of 
these  methods  is  given  in  [11]  .  Note  that  in  all  cases, 
conservativeness  of  the  schema  is  an  ever-present  requirement. 
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